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Under laboratory conditions coronaviruses were shown to have a high frequency of recombination. 
In The Netherlands, vaccination against infectious bronchitis virus (IBV) is performed with vaccines 
that contain several life-attenuated virus strains. These highly effective vaccines may create ideal 
conditions for recombination, and could therefore be dangerous in the long term. This paper 
addresses the question of the frequency of recombination of avian coronavirus IBV in the field. A 
method was sought to detect and quantify recombination from sequence data. Nucleotide sequences 
of eight IBV isolates in a region of the genome suspected to contain recombination, were aligned 
and compared. Phylogenetic trees were constructed for different sections of this region. Differences 
in topology between these trees were observed, suggesting that in three out of eight strains in vivo 


RNA recombinant had occurred. 
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Introduction 


Infectious bronchitis virus (IBV) causes considerable 
damage in the poultry industry by infections of the 
respiratory tract and the reproductive organs’. IBV is 
the prototype strain of the Coronaviridae, a family of 
positive stranded RNA viruses with a genome of about 
28kb. This genome contains the information for three 
structural proteins and a number of enzymes involved in 
virus replication (Figure 2). Although vaccines offer 
protection, they are often made ineffective by the 
continuous emergence of new serotypes?*. Serotype- 
specific protective immunity is thought to be mediated 
by antibodies to a protein on the surface of the virion, 
the so-called peplomer protein*~’. The peplomer protein 
is synthesized as a precursor glycoprotein, which is 
cleaved into the subunits Si and S2, derived from its 
N- and C-terminal half, respectively*~*°. Point mutations 
are obviously involved in the generation of new serotypes. 
However, previous studies’! suggested that RNA recom- 
bination is also involved in the generation of new 
antigenic variants. To address the question of the 
frequency of recombination in the generation of new field 
isolates the nucleotide sequences in five windows of 
homologous sequences of the genomes of eight IBV 
strains have been compared. To make this comparison 
the sequences of the S2 genes of D207 and D1466 were 
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newly determined. Both strains are Dutch field isolates 
and are included in vaccines currently used in The 
Netherlands. The alignment of the sequences suggests 
that the Dutch field strain D207!?, the British field strain 
6/823 and the Japanese field strain KB8523'* result from 
RNA recombination events. 


Materials and methods 


The origin of the IBV strains as well as the procedures 
used for the cloning and sequencing of the peplomer 
genes have been published previously***:1°. Sequence 
data were processed using the programs of Microgenic 
(Beckman Instruments, Ref. 16). Phylogenetic trees were 
derived from nucleotide difference-matrices using the 
program of Dr J. Felsenstein (University of Seattle) 
distributed as part of the PHYLIP package 2.6'’. 


Results and discussion 


The S2 sequences 


The newly determined D207 and D1466 sequences are 
based on data from two or more independent cDNA 
clones. With strain D1466, two nucleotide differences 
between cDNA clones were observed. Neither of these 
differences resulted in an amino acid replacement. From 
the nucleotide sequences (submitted to the EMBL 
Genbank and DDBJ Nucleotide Sequence Databases), 
amino acid sequences have been deduced. In Figure 1 
the amino acid sequences are listed together with those 
of all other known S2 sequences. 

Assuming the conserved Arg-Arg-X-Arg-Arg-Ser se- 
quence to be the cleavage site between S1 and S2%, the 
Ser at position 546 (Figure 1) will become the N-terminal 
S2 residue. In contrast to S1, where many insertions/ 
deletions were found, the variation in length of S2 is 
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Figure 1 Amino acid sequences of IBV S2 proteins. The sequence of M41 variant M41U has been listed completely; from other strains only the 
differences with this sequence are shown. The sequences of M4iU, and of the M42 variant M42S are from Kusters et a/."". The sequence of M42 
variant M42H is from Binns et a/.?°. The sequences of M41 variant M41H and of strain 6/82 are from Binns et a/.*. The sequence of KB8523 is from 
Sutou et a/.*. Potential N-glycosylation sites (NXS or NXT, except when X=P) are underlined. Asterisks indicate the stop codons 


caused solely by extension of the open reading frame. 
Substitution of the T in the TAA stop codon of strains 
M41H and D1466 by a G in all other strains resulted in 
an extension of the S2 protein by nine amino acid 
residues. Generally the amino acid sequences of the S2 
proteins differ by less than 10%, strain D1466 being the 
exception (Table 1). However, 19 out of 20 cysteine 
residues and most glycosylation sites are completely 
conserved. 


Analysis of nucleotide sequences for recombination 


The nucleotide and amino acid sequences were 
compared for possible recombination points. From this 
comparison we conclude that the Japanese field strain 
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KB8523 is probably a recombinant. The S2 protein of 
this strain is almost identical to S2 of M41, except for 
the region corresponding to amino acids 929 to 1102 
(Figure 1) where the sequence is remarkably similar to 
D1466. This suggests a recombination event with two 
crossovers in this region. The alignment of the S2 genes 
indicates a second potential recombinant: the British field 
isolate 6/82. Its S1 gene and most of the S2 gene are 
almost identical to those of strain D207. At the 3’ end 
of the S2 sequence (at the Lys-1105 codon) the percentage 
of identical nucleotides switches from 99% upstream to 
57% downstream. At exactly the same position the 
nucleotide sequence of 6/82 becomes highly similar to 
the D1466 sequence: from 73% identity upstream to 97% 
downstream. The nucleotide identity between D207 and 
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Table 1 Difference matrix of S2 sequences of IBV strains* 


Virus M41U M41H M42H M42S KB8523 6/82 D207 D1466 
i a gn ce a eee Ce 
M41U - 1.4 2.9 2.6 6.7 7.0 6.9 26.7 
M41H 1.4 - 4.2 3.8 8.0 8.3 8.2 25.6 
M42H 2.9 4.2 = V1 7.8 7.5 6.9 27.4 
M42S 2.6 3.8 11 - 7.5 7.5 7.2 27.4 
KB8523 6.7 8.0 7.8 7.5 a 9.0 8.8 26.6 
6/82 7.0 8.3 7.5 7.5 9.0 - 1.8 25.4 
D207 6.9 8.2 6.9 7.2 8.8 1.8 = 26.6 
D1466 26.7 25.6 27.4 27.4 26.6 25.4 26.6 - 


*The figures represent the percentages of non-identical amino acids. Designation of strains is as in Figure 1 
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Figure 2. Phylogenetic relationships of IBV strains listed in Figure 1 (A) Schematic representation of the genome of IBV. The roman numerals |-V 
indicate the regions used to construct the phylogenetic trees. (B) Phylogenetic trees for the regions I-V as indicated in (A). |, Tree of S1 sequences; 
Il, tree of S2 sequences upstream residue 929 (Figure 1); Ill, tree of S2 sequences between residues 929 and 1105; IV, tree of S2 sequences downstream 
residue 1105 until the stop codon on position 1182; V, tree of sequences from the 5’ part of the E1 gene and part of the preceding intergenic sequences”’. The 
figures indicate the permillage of replaced nucleotides. The trees are unrooted, i.e. the position of the hypothetical ancestral IBV sequences are 
unknown. For the construction of the trees of this figure, only the M41H and M42U sequences were used. Replacing them by M41U and/or M42S, 
respectively, did not significantly alter the topology of the trees (data not shown) 


D1466 is about the same on both sides of this putative these results more rigorously. This was achieved by 
crossover site: 73 and 67%, respectively. comparing phylogenetic trees for each of five sections 

The above comparisons are both time consuming and (windows) of the S2 gene and its flanking sequences (see 
subjective. A method was sought to test and quantify Figure 2a). The S2 gene was divided into three sections. 
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The first section (II in Figure 2b) contained the 5’-part 
of the S2 gene up to the putative KB8523 crossover site 
(in the codon for Ile-929 in Figure 1). A second S2 tree 
(III in Figure 2b) was constructed for the part of the S2 
gene between the KB8523 and 6/82 crossover sites (from 
Ile-929 to Lys-1105). The 3’-border of this region did not 
exactly coincide with the second KB8523 crossover site 
(in the Ser-1102 instead of the Lys-1105 codon) but 
changing this border had no significant effect on the 
topology. The third S2 tree (IV in Figure 2b) was 
calculated for the region from the 6/82 crossover site (in 
the Lys-1105 codon) until the stop codon on position 
1182. 

Sequences upstream and downstream from the S2 gene 
were also used to construct phylogenetic trees. The tree 
based on S1 gene sequences (I in Figure 2b) has the same 
topology as an S2 tree based on region II in S2. KB8523 
has shifted to the D1466 branch in a tree for region III of 
S2. Region IV of S2 yields a tree in which 6/82 has shifted 
from the D207 to the D1466 branch, while KB8523 is 
again near M41 and M42. Finally a tree (V) constructed 
from the El gene sequences shows that, compared to tree 
IV, also strain D207 has shifted towards the D1466 
branch, and therefore may also be a recombinant. 


Conclusions 


With the murine coronavirus MHV a high frequency of 
recombination was found both in vitro and in mouse 
brain after infection with a ts mutant of MHV strain A59 
and wild type JHM-virus!®!°. 

Phylogenetic trees constructed from five different 
windows of the genomes of eight IBV strains were used 
to detect crossover events. A clear advantage of this 
method is that, due to its simplicity, large numbers of 
homologous genes can easily be screened. The data 
presented in this paper suggest that genomic recombi- 
nation not only occurs under laboratory conditions but 
also plays an important role in the generation of new 
virulent strains; at least three IBV strains, all isolated 
from outbreaks, are potential recombinants. Previous 
data'’ already indicated that the Dutch field isolate 
V1397 is also the product of a recombination. Considering 
the limited number of sequences analysed so far, in vivo 
recombination of IBV seems to occur at a rather high 
frequency. This has implications for the use of IBV 
vaccines consisting of several serotypes of live attenuated 
IBV strains. These vaccines offer the best protection, but 
they may prove dangerous in the long term by the 
induction of new variants, not only by point mutations 
but also through recombination. In this respect vacci- 
nation protocols that combine IBV with other live 
attenuated viruses are equally suspect. They may create 
ideal conditions for non-homologous recombination 
which may result in the emergence of recombinant virus 
species with unwanted properties. 
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